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Abstract 

Background: The single nucleotide polymorphism (SNP) rs2615977 is associated with osteoarthritis (OA) and is 
located in intron 31 of C0L1 1A1, a strong candidate gene for this degenerative musculoskeletal disease. 
Furthermore, the common non-synonymous C0L1 1A1 SNP rsl 676486 is associated with another degenerative 
musculoskeletal disease, lumbar disc herniation (LDH). rsl 676486 is a C-T transition mediating its affect on LDH 
susceptibility by modulating C0L11A1 expression. The risk T-al lele of rsl 676486 leads to reduced expression of the 
COLII A) transcript, a phenomenon known as allelic expression imbalance (AEI). We were keen therefore to assess 
whether the effect that rsl 676486 has on C0L1 1A1 expression in LDH is also observed in OA and whether the 
rs261 5977 association to OA also marked AEI. 

Methods: Using RNA from OA cartilage, we assessed whether either SNP correlated with C0L11A1 AEI by 1) measuring 
C0L11A1 expression and stratifying the data by genotype at each SNP; and 2) quantifying the mRNA transcribed 
from each allele of the two SNPs. We also assessed whether rsl 676486 was associated with OA susceptibility using 
a case-control cohort of over 18,000 individuals. 

Results: We observed significant AEI at rsl 676486 (p< 0.0001) with the T-allele correlating with reduced COL11A1 
expression. This corresponded with observations in LDH but the SNP was not associated with OA. We did not 
observe AEI at rs2615977. 

Conclusions: COL11A1 is subject to AEI in OA cartilage. AEI at rsl 676486 is a risk factor for LDH, but not for OA. 

These two diseases therefore share a common functional phenotype, namely AEI of COL11A1, but this appears to 
be a disease risk only in LDH. Other functional effects on COL11A1 presumably account for the OA susceptibility 
that maps to this gene. 
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Background 

Primary osteoarthritis (OA) is a common disease of 
the synovial joint characterised by articular cartilage 
thinning and loss, which is often accompanied by al- 
teration in the normal function of other tissues of 
the joint [1]. A number of epidemiological studies 
have demonstrated that osteoarthritis has a major 
genetic component, although elucidation of this has 
proven challenging [2]. 
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Genome-wide association scans (GWASs) are now 
being employed as a powerful, objective tool for map- 
ping OA susceptibility loci [3-6]. In some circum- 
stances the resulting signals are to regions of the 
genome that do not harbour obvious OA candidates, 
whereas in others the signals are close to, or even 
within, a gene that would be considered a strong can- 
didate. One example of the latter emerged from the 
original arcOGEN GWAS [4]. This GWAS and its 
subsequent replication analysis was performed on in- 
dividuals of European descent and included the full 
arcOGEN samples and, amongst other signals, identi- 
fied association to single nucleotide polymorphism 
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(SNP) rs2615977, with a p-value of 1.2xl0~ 5 and an 
odds ratio of 1.10 for the T- allele of the SNP. 

rs2615977 resides within intron 31 of COL11A1. 
This gene codes for the al polypeptide chain of type 
XI collagen, a quantitatively minor but structurally 
important component of the cartilage extracellular 
matrix. Mutations within COL11A1, typically involv- 
ing the substitution of highly conserved glycine resi- 
dues or the deletion of whole exons, are responsible 
for certain forms of Stickler syndrome, a rare and 
usually dominant Mendelian condition that presents 
with a severe and premature form of OA manifesting 
during adolescence [7,8]. 

Lumbar disc herniation (LDH) and lumbar disc de- 
generation (LDD) can be viewed as having similar 
aetiological routes to OA, with each condition being 
age-associated and involving the degeneration of cartil- 
age, albeit of different types, with OA being hyaline 
and LDH/LDD being fibrocartilage. The common non- 
synonymous SNP rs 1676486, which is located in exon 
62 of COL11A1, has been reported to associate with 
LDH in a Japanese population [9]. rs 1676486 is a C-T 
transition that results in a proline to serine substitu- 
tion. Using an in vitro assay the investigators were able 
to demonstrate that the LDH-associated T-allele of 
rsl676486 correlated with decreased stability of the 
COL11A1 transcript relative to the C-allele. Such a differ- 
ence in allelic expression, whether mediated by differen- 
tial mRNA transcription or mRNA transcript stability, is 
known as allelic expression imbalance (AEI). As such, 
the investigators concluded that a quantitative deficiency 
of COL11A1, and by extrapolation a quantitative defi- 
ciency of the al polypeptide chain of type XI collagen, 
was a contributor to LDH susceptibility. 

As more susceptibility loci are identified for common 
human diseases it is becoming clear that many, if not the 
majority, of associated alleles contribute to disease risk 
by influencing gene expression, typically by modulating 
rates of transcription or transcript stability [10,11]. The 
T-allele of rsl676486 is an example of such an effect in 
LDH. It is also becoming apparent that in most of these 
cases the effect that the susceptibility allele has on the tran- 
script can be organ, tissue or even cell type-specific [12,13]. 
It is important therefore to study RNA extracted from tis- 
sues relevant to the disease of interest when assessing the 
potential effect of a functional susceptibility allele. 

Based on these observations we hypothesised that the 
OA association to COL11A1, marked by rs2615977, may 
be mediated by modulating the expression of the gene in 
articular cartilage. In addition, we hypothesised that the 
effect that rsl676486 has on COL11A1 expression in 
LDH may also be observed in OA. To investigate these 
hypotheses we have quantitatively measured overall ex- 
pression of COL11A1 in cartilage and then stratified our 



data by genotype at rs2615977 and at rsl676486. We 
have also tested for AEI of COL11A1 using assays that 
can accurately discriminate and quantify the mRNA out- 
put from each allele of a transcript SNP. 

Methods 

Patients and tissues 

Macroscopically normal articular cartilage tissue located 
away from the OA lesion was obtained from individuals 
undergoing elective joint replacement for OA of the hip 
(total hip replacement, THR) or of the knee (total knee 
replacement, TKR), as described in detail previously 
[14,15]. The Newcastle and North Tyneside research eth- 
ics committee granted ethical approval for the collection 
(REC reference number 09/H0906/72) and written in- 
formed consent was obtained from donors for the use of 
their tissue, and permission for publication of their age 
and sex. Details regarding the patients can be found in 
Additional file 1: Table SI and Additional file 2: Table S2. 

Nucleic acid extraction 

On the day of surgery, the tissue specimens were snap- 
frozen at -80°C. For each individual tissue sample, 0.5- 
1.0 g of frozen tissue was ground to a powder using a 
Retsch mixermill 200 (Retsch Limited, Leeds, UK) under 
liquid nitrogen, which causes the sample to become brit- 
tle and prevents the RNA from degrading. Genomic 
DNA and RNA were then extracted from the ground 
tissue samples using an EZNA DNA/RNA Isolation Kit 
and a protocol established for human tissue (Omega 
bio-tek, R6731-02). The nucleic acid was quantified using 
a NanoDrop ND-1000 Spectrophotometer (NanoDrop 
Technologies, Wilmington, USA). 

SNPs 

We studied three COL11A1 SNPs: the OA associated 
SNP rs2615977, which is located in intron 31; the LDH 
associated SNP rsl676486, which is located in exon 62; 
and SNP rs9659030, which is located in the 3'UTR 
(Table 1). rs9659030 was analysed for AEI in the absence 
of a COL11A1 transcript SNP in high linkage disequilib- 
rium (LD) with the OA SNP rs2615977. rs2615977 is 
98.2 kb from rsl676486 and 110 kb from rs9659030. 

SNP genotyping 

The three COL11A1 SNPs were genotyped by restriction 
fragment length polymorphism (RFLP) analysis (patient 
genotypes are listed in Additional file 1: Table SI). The 
primer sequences and restriction enzymes used can be 
found in Additional file 3: Table S3. 

cDNA synthesis and quantitative real-time PCR 

RNA extracted from cartilage was reverse transcribed 
using the Superscript First-Strand cDNA synthesis kit 
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Table 1 The three C0L11A1 SNPs studied in this report 

SNP Alleles MAF 1 Pairwise linkage disequilibrium 

Relative to rs2615977 Relative to rs1 676486 Relative to rs9659030 

ry ? ry ? ry ^ 

rs2615977 T/G 0.239 - - 0.58 0.03 0.39 0.08 

rsl 676486 C/T 0.202 0.58 0.03 0.33 0.09 

rs9659030 A/G 0.173 0.39 0.08 0.33 0.09 
1 Minor allele frequency in Europeans (HapMap CEU). 



(Invitrogen). Gene expression was measured by quan- 
titative real-time PCR using PrimeTime Mini qPCR 
Assays (Integrated DNA Technologies, Iowa, USA) 
and an ABI PRISM 7900HT Sequence Detection Sys- 
tem (Applied Biosystems). COL11A1 expression was 
measured relative to the average expression of the 
housekeeping genes HPRT1, GAPDH and 18S: for 
each cDNA sample three pipetting replicates we 
performed per gene and the nine housekeeper cycle 
threshold (Ct) values (3x3) were combined to derive 
the mean Ct which was used as the control to compare 
against COL11A1. Expression of COL11A1 relative to 
the housekeepers was determined using the 2~ A ct 
method. Mann- Whitney U and Kruskal-Wallis tests 
were performed to assess whether COL11A1 expression 
relative to genotype at rs2615977 and genotype at 
rsl676486 differed significantly from the null. All pri- 
mer and probe sequences can be found in Additional 
file 4: Table S4. 

AEI analysis 

AEI was assessed using the transcript SNPs rsl676486 
and rs9659030 and quantitative real time PCR geno- 
typing assays, purchased from Applied Biosystems. 
These assays are standard real time assays except that 
they employ a probe (labeled with FAM or VIC) spe- 
cific to each of the two alleles of a SNP. Real time 
PCR was carried out according to the manufacturer's 
instructions. In brief, 10 ng of cDNA or 20 ng of 
DNA extracted from the cartilage was added to 5 ul 
of TaqMan Universal Master Mix II no UNG (Applied 
Biosystems) and 0.25 ul of 40 x TaqMan assay in a 
10 ul reaction. The samples were then amplified on 
an ABI PRISM 7900HT Sequence Detection System. 
The reactions were performed with five pipetting rep- 
licates and the allelic ratios were calculated using the 
formula (2' FAM Ct )/(2" VIC ct ). 

For each assay, the ratios between the amounts of 
each allele in every sample were calculated for genomic 
DNA and cDNA. For each sample the average allelic 
ratio for genomic DNA, which represents the 1:1 ratio, 
was then used to normalise the cDNA ratio to generate 
a corrected allelic ratio. To determine if there was an 
overall difference in expression between alleles for a 



particular tissue across all patients the mean allelic ra- 
tios for the patient cDNAs were compared to the mean 
allelic ratios for the patient genomic DNAs using a 
Mann-Whitney U test. 

Quality control checks on the AEI assays 

The capacity for each assay to discriminate between 
SNP alleles was verified using a standard curve of artifi- 
cially created allelic ratios using DNA from samples 
homozygous for the major and the minor alleles. To as- 
sess the amplification efficiency of each AEI probe (two 
per assay) a standard curve was produced by performing 
AEI on a serial dilution of 200 ng to 1.5635 ng of DNA 
from a heterozygous sample. 

Association of rs 1676486 with OA 

We analysed the LDH associated SNP rsl676486 for evi- 
dence of association to OA using data from the UK 
arcOGEN study. This is a GWAS performed for 485,491 
autosomal SNPs and conducted on 7,410 cases with se- 
vere OA of the hip or the knee, 80% of whom had 
undergone total joint replacement, and on 11,009 publi- 
cally available population controls [6]. 

Results 

Quantitative expression of COL11A1 in cartilage stratified 
by genotype at the OA and LDH associated SNPs 

To initially assess the degree to which rs2615977 and 
rsl676486 correlate with expression of COL11A1, we 
used quantitative real-time PCR to measure the level 
of expression of the gene in cDNA that had been 
synthesised from cartilage RNA. We then stratified 
the data by genotype at rs2615977 and by genotype 
at rsl676486 (Figure 1). Only one patient was homo- 
zygous for the minor allele of rsl676486. This is a 
reflection of the low frequency of TT homozygotes 
for this SNP, with a frequency of only 1.8% reported 
in HapMap CEU. As such, for this SNP a compari- 
son was made between those patients who carry the 
LDH risk T-allele and those who do not, that is, be- 
tween CC homozygotes versus CT heterozygotes and 
TT homozygotes combined. There was no significant 
correlation (p < 0.05) between the level of COL11A1 
expression and the genotype at either rs2615977 or 
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Figure 1 Columnar scatter plot of the quantitative expression of COL11A1 in cartilage cDNA stratified by genotype at A) the OA 
associated SNP rs2615977 and B) the LDH associated SNP rs1 676486. Due to the low frequency of TT homozygotes at rs1 676486 the 
analysis in B) was between CC homozygotes and T-allele carriers (CT and TT combined) at this SNP. n is the number of patients studied for each 
group. The horizontal lines in each plot represent the mean and the standard error of the mean. P-values were calculated using a Kruskal-Wallis 
test for A) and a Mann-Whitney U test for B). 



rsl676486. We subsequently stratified the patients 
into those who had undergone TKR and those who 
had undergone THR, and repeated the analysis. Again, 
there was no correlation between COL11A1 expression 
and genotype at either SNP (Figure 2). 



AEI analysis 

The above analysis did not provide evidence for a correl- 
ation between COL11A1 expression and genotype at 
rs2615977 or genotype at rsl676486. However, this tech- 
nique is vulnerable to the natural fluctuation in gene 
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Figure 2 Columnar scatter plots of the quantitative expression of COL11A1 in OA cartilage cDNA stratified by A) genotype at the LDH 
associated SNP rs1 676486 in knee cases, B) genotype at the LDH associated SNP rs1 676486 in hip cases, C) genotype at the OA 
associated SNP rs2615977 in knee cases, and D) genotype at the OA associated SNP rs2615977 in hip cases. Due to the low frequency of 
TT homozygotes at rs 1 676486 the analysis in A) and B) was between CC homozygotes and T-allele carriers (CT and TT combined) at this SNP. 
Due to the absence in the hip strata of GG homozygotes at rs261 5977 the analysis in D) was between TT homozygotes and TG heterozygotes at 
this SNP. n is the number of patients studied for each group. The horizontal lines in each plot represent the mean and the standard error of the 
mean. P-values were calculated using a Mann-Whitney U test for A), B) and D) and a Kruskal-Wallis test for C). 
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expression, which can decrease its sensitivity and ac- 
curacy [16]. For example, it is apparent from 
Figures 1 and 2 that between individuals of the same 
genotype the expression of C0L11A1 can vary by 
greater than 200 fold. As such, quantitative PCR can 
generate false negative results when used in this 
manner. An alternative method for correlating geno- 
type with expression levels that is not vulnerable to 
such fluctuation is allelic expression imbalance, ab- 
breviated to AEI and also known as differential allelic 
expression. This entails using a transcript SNP to dif- 
ferentiate and measure the cDNA output from the 
two alleles of a heterozygote. Genomic DNA is used 
as the comparator control to provide the 1:1 ratio to 
which the cDNA allelic ratio is compared. We there- 
fore assessed whether COL11A1 demonstrated AEI 
and if so, whether this correlated with genotype at 
rs2615977 or at rsl676486. 

rsl676486 is a transcript SNP and it could itself be 
used in the AEI analysis. However, rs2615977 is intronic 
and could not be used. A search of the Ensembl data- 
base (http://www.ensembl.org/) for COL11A1 transcript 
SNPs in high LD (r 2 > 0.5) with rs2615977 identified only 
two SNPs, which happen to be in perfect LD (r 2 = 1, 
D'=l) with the SNP: rsl012281 and rsl012282, both 
within exon 45. However, according to the UCSC data- 
base (http://genome.ucsc.edu/), the sequence harbouring 
these SNPs is not part of the COL11A1 transcript. To re- 
solve this discrepancy we amplified across the region 
using OA cartilage cDNAs and genomic DNAs as tem- 
plates. The same primers were used for the two tem- 
plates and we successfully amplified DNA but we failed 
to detect a cDNA product (results not shown). We con- 
cluded therefore that the UCSC database was accurate 
and that rsl012281 and rsl012282 are not transcript 
SNPs and could not therefore be used. Since there are 
no transcript SNPs in high LD with rs2615977 we opted 
to use a transcript SNP that had a high heterozygosity 
value, thus ensuring that a large number of our patients 
would be informative for the AEI analysis if we focused 
on individuals who were compound heterozygotes for 
rs2615977 and the transcript SNP. Due to the small size 
of the majority of COL11A1 exons we needed to avoid a 
transcript SNP that was close to an intron/exon bound- 
ary, as this would not allow us to use the same primers 
and probe when targeting cDNA and DNA in our AEI 
analysis. The transcript SNP that we chose was the 
3'UTR SNP rs9659030, which has a heterozygosity value 
of 0.31 in Europeans and is over 350 bp from the nearest 
exon boundary. Although rs9659030 is not in LD with 
rs2615977 (pair wise r 2 = 0.08, pair wise D' = 0.40), by 
studying compound heterozygotes we were able to use 
rs9659030 to functionally separate the two rs2615977 al- 
leles. Although we could not determine phase between 



rs2615977 and rs9659030 we would expect a high pro- 
portion of the compound heterozygotes to demonstrate 
AEI if rs26 15977 is itself, or if it correlates with, a 
polymorphism that influences gene expression at an al- 
lelic level. 

rsl 676486 

We carried out AEI analysis on 22 OA patients hetero- 
zygous for the LDH SNP rsl676486, all of whom 
showed a significant (p < 0.05) decrease in expression of 
the T-allele relative to the C-allele with an average al- 
lelic ratio of 0.36 (p < 0.0001) (Figure 3 and Additional 
file 5: Table S5). That is, for every one C-allele cDNA 
molecule there were only 0.36 cDNA molecules of a 
T-allele. It was the T-allele of rsl676486 that showed 
reduced stability in the original LDH study [9]. Our 
functional data therefore supported the functional data 
from the LDH study. 

rs9659030 

We carried out AEI analysis on 26 OA patients hetero- 
zygous for rs9659030 (Figure 4 and Additional file 5: 
Table S5). Thirteen of the 26 were compound heterozy- 
gous for rs2615977 and rs9659030 but only one of the 
13, patient 20, showed a significant AEI. The two alleles 
of the OA associated SNP rs2615977 do not therefore 
correlate with differences in expression of COL11A1. In 
addition to patient 20 three other patients also showed a 
significant AEI (patients 13, 37 and 39 in Figure 4). 
These three patients were not heterozygous at 
rs2615977 but they were heterozygous for the LDH SNP 
rsl676486 and it is presumably this that accounts for 
the AEI observed in these patients. 

We were intrigued to note that the further nine 
patients who were also compound heterozygous for 
rsl676486 and rs9659030 did not demonstrate sig- 
nificant AEI when their allelic expression was mea- 
sured using rs9659030 (patients 1, 3, 8, 10, 11 24, 
31, 46 and 76 in Figure 4) - based on the AEI re- 
sults for rsl676486 reported earlier we had antici- 
pated that all nine would have demonstrated 
significant AEI with rs9659030, particularly in light 
of the fact that eight of the nine had been tested 
with rsl676486 and had shown significant AEI with 
this SNP (patients 1, 3, 8, 10, 11, 24, 31 and 46 in 
Figure 3). Quality control (QC) checks on the 
rs9659030 AEI assay revealed that it is able to accur- 
ately discriminate between and measure differences 
in the levels of the two alleles of this SNP (Add- 
itional file 6: Figure SI) and that this capacity is well 
within the range of AEI observed at rsl676486. For 
example, the average rsl676486 C/T ratio for the 22 
patients of Figure 3 and Additional file 5: Table S5 
was 0.36 and it is clear from Additional file 6: Figure 
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Figure 3 Allelic expression analysis in cartilage samples from OA patients assessed using the LDH associated SNP rsl 676486. The 

genotypes of the patients at rs1676486, at the OA associated SNP rs261 5977 and at the 3'UTR SNP rs9659030 are shown. P-values were calculated 
by comparing the cDNA allelic ratio to the DNA allelic ratio using a Mann-Whitney U test. ** denotes p < 0.01. 



S1A and SIB that the rs9659030 AEI assay would 
have detected such a difference if it existed at this 
SNP. These QC checks confirm that the rs9659030 
AEI assay is robust. Overall, these results imply that 
there is variation in the degradation rate of the 
COL11A1 transcript across its length, which would 
mean that the AEI results obtained could partly be 
dependent on the physical position in the transcript 
of the SNP used to measure AEI. 



Genetic association analysis of rsl 676486 with OA 

We finally assessed whether rsl676486 was associated to 
OA using the arcOGEN GWAS dataset. This revealed 
that the SNP showed only marginal evidence of associ- 
ation, with a p-value of 0.012 in the comparison of 7,410 
cases and 11,009 controls (Additional file 7: Table S6). 
The association signal was in the same direction as seen 
for the Japanese LDH study, in that it was the T- allele of 
rsl676486 that was elevated in cases (frequency of 0.201) 
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Figure 4 Allelic expression analysis in cartilage samples from OA patients assessed using the 3'UTR SNP rs9659030. The genotypes of 
the patients at rs9659030, at the OA associated SNP rs261 5977 and at the LDH associated SNP rs1676486 are shown. P-values were calculated by 
comparing the cDNA allelic ratio to the DNA allelic ratio using a Mann-Whitney U test. ** denotes p < 0.01 . 



Raine et al. BMC Musculoskeletal Disorders 2013, 14:85 
http://www.biomedcentral.com/1471-2474/14/85 



Page 7 of 8 



versus controls (frequency of 0.1905) with an odds ratio 
of 1.069 (95% confidence interval of 1.015-1.127). Strati- 
fication of the arcOGEN data by site of OA (hip or knee) 
did not enhance the association signal (Additional file 7: 
Table S6) nor was there a more significant signal when 
the data was stratified by sex, or by sex combined with 
site (data not shown). 

Discussion 

Allelic expression imbalance (AEI) has been shown to be 
an important mechanism through which OA susceptibil- 
ity acts, with an excellent example being the SNP 
rsl43383 located in the 5'UTR of the growth differenti- 
ation factor 5 gene GDF5; this SNP is robustly associated 
with knee OA across various ethnic groups, with the 
OA-associated T-allele of the polymorphism mediating 
reduced expression of GDF5 in cartilage and other joint 
tissues [reviewed in 2]. Based on observations such as this, 
we were keen to assess whether the COL11A1 OA associ- 
ation, marked by rs2615977, also correlated with AEI in 
OA cartilage. A previously published study of LDH dem- 
onstrated that this gene was subject to AEI and we there- 
fore also assessed whether the LDH SNP rsl676486 
correlated with AEI in OA cartilage and if so, whether the 
polymorphism was genetically associated with the disease. 

Our results did not provide evidence for a correlation 
between the OA associated SNP rs26 15977 and 
COL11A1 AEI but they did provide compelling evidence 
for a correlation between the LDH associated SNP 
rsl676486 and AEI at this gene. As for LDH, it was the 
T-allele of rsl676486 that correlated with reduced ex- 
pression of COL11A1. However, the polymorphism dem- 
onstrated only very weak evidence of association with 
OA and as this association data derives from a GWAS of 
over 480,000 SNPs, the very moderate p-value does not, 
in our view, support rsl676486 as an OA susceptibility 
locus despite its functional effect. From our rsl676486 
data we therefore conclude that in the cartilage of OA 
patients who have undergone elective hip or knee joint 
replacement surgery, COL11A1 is subject to AEI, but that 
this is not a risk factor for OA. Mature hyaline cartilage 
tissue is therefore resilient to such COL11A1 allelic im- 
balance effects, unlike the fibrocartilage of LDH, which is 
vulnerable to COL11A1 AEI. 

The lack of a correlation between the OA associated 
SNP rs2615977 and COL11A1 expression implies that the 
association marked by this SNP is operating by a route 
other than an effect on the expression of the gene in ma- 
ture cartilage. However, it is possible that we failed to de- 
tect any such correlation, since in our AEI analysis we used 
a SNP, rs9659030, that is not physically close to rs2615977. 
After all, our comparison of AEI in individuals compound 
heterozygous for rs9659030 and rsl676486 did imply that 
AEI may not be uniform across the COL11A1 transcript. 



However, the fact that AEI in mature cartilage as mea- 
sured with rsl676486 is not an OA risk factor makes it 
seem unlikely that the OA risk marked by rs2615977 
would be operating via this mechanism in this tissue. It 
is possible that rs2615977, or a SNP in high LD with it, 
is modulating the expression of COL11A1 at time points 
not studied here, for example during skeletogenesis or in 
adolescent tissues. Of possible relevance to this is the 
presence in the ENCODE database [17] of SNP 
rsl 1164648, which is in perfect LD with rs2615977 and 
which is located within a CTCF-binding site that is pre- 
dicted to act as a transcriptional insulator. CTCF is a 
transcriptional repressor and CTCF-dependent insulators 
have been shown to regulate gene expression by modulat- 
ing chromatin conformation during development [18]. If 
cartilage tissues were to become available from an ad- 
equate number of young donors then assessing COL11A1 
AEI relative to genotype at rs2615977 or rsl 1164648 may 
be insightful. 

As noted in the introduction, Stickler syndrome is a 
rare and usually dominant condition presenting with a 
severe and early onset form of OA [7,8]. Stickler syn- 
drome can present with a high degree of inter and 
intra-familial variability and as such it may be of bene- 
fit to assess whether the possession of the low express- 
ing T-allele of rsl676486, which we have shown is 
functionally active in cartilage and which has a fre- 
quency of >20% in Europeans and Asians, correlates in 
any way with this phenotypic heterogeneity; possessing 
a low expressing COL11A1 allele alongside a mutant 
allele may exacerbate the phenotype whilst if the mu- 
tant COL11A1 allele is on a low- expressing allele then 
this could attenuate the phenotype. 

Conclusions 

Overall our study has identified that COL11A1 AEI is 
common to both LDH and OA, but that it is not a risk 
factor for OA. Presumably other functional effects must 
be accounting for the OA associated signal at this gene. 
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